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We consider the two dimensional disordered Bose gas which present a metallic state at low tem- 
peratures. A simple model of an interacting Bose system in a random field is propose to consider 
the interaction effect on the transition in the metallic state. 



I. INTRODUCTION 

The recent experimental data obtained on the two- 
dimensional (2d) systems [H, [2, [l] showed the possibility 
of the occurrence for the metallic phase in the insulator- 
superconductor transition. The model proposed by Das 
and Doniach [3| is a uniform phase lacking both phase 
and charge order. This phase presents translationally 
and rotationally invariance. It has to appear at finite 
temperature, because at T = there is no phase tran- 
sition which is lacking both charge and phase order 
The occurrence of the metallic phase glass in a disordered 
two dimensional {2d) bosonic system has been predicted 
by Wagenblastl et al. Q in the insulator-superconductor 
transition, but only at the separatrix line. The important 
result has been obtained by Dalidovich and Phillips Q, 
who shoved that this state can appear in disordered two 
dimensional {2d) Bose systems. This model is based on 
the quantum version of the spin glass state which can be 
described using the replica trick method. The initial free 
energy has the universality class z — 1 {z is the dynami- 
cal critical exponent) with a local breaking spin rotation 
invariance. The effective order parameter , is given in 
terms of the phase (p{r, r) (t is the imaginary time) as 
, Q{t) =< exp(i(p(ry , 0) exp(«(/j(rjj, t))) > and relaxes 
exponentially fast to its equilibrium value. The quantum 
decay of Q{t) is like and this gives rise to excitations 
which scales as and changes the universality class to 
z = 2 [1]. Such a dissipative model has been treated in 
in the Gaussian approximation and the conductivity 
a{T = 0, w = 0) has been obtained as finite, which proves 
that the state is metallic. 

In this paper we present an equivalent model for the 
2d disordered Bose system which can be treated with the 
Renormalization Group (RG) method. The model was 
used to describe the spin glass state where disor- 

der was introduced by a random field which can be con- 
nected to the effective order parameter of the spin glas 
state. The quantum version of this method was given in 
[ill [T^ where was applied for the study of the critical 
behavior of the interacting Bose system. However, we 
mention that the theory of the critical behavior of the 



transport is not trivial because at finite temperature the 
limits a{T = 0, uj) and a{u) = 0, T) are not equal [Hill, 
but we will be interested in the calculation of the ct(0, 0) 
near the critical region. The model of our disordered 
Bose system will be presented in Sec. II, where we will 
discuss the equivalence with the model from Ref . [§| . In 
Sec. Ill we solve the flow equations of the RG. equations 
of this model. The strategy of our calculations is dif- 
ferent to this from [111, [13] , because we are interested in 
the case of d = 2 and z = 2. In this case we can solve 
the differential equations exactly, studding the influence 
of the temperature on the transition. In fact this pro- 
cedure give us the possibility to analyze the relevance of 
the parameters T and the interaction uq near the fixed 
point T = 0. We expect that the new result for conduc- 
tivity to contain a correction given by the interactions 
between bosons which is in fact small and are in fact 
corrections to the Gaussian behavior. The conductivity 
near the critical region will be calculated following the 
method given in 0, and applied recently [lB| for the 
study of the study of the fluctuations of the conductivity 
in the d-wave superconductors near the critical disorder. 
The last section. Sec. IV, will be devoted to the calcula- 
tion of the conductivity and the result will be compared 
with the simple Gaussian model. . The values given by 
the experimental are different from the standard result 
obtained for the simple model which gives ct(0) = 4e^//i 
( 2e is the pair charge and h is the Planck constant) 
the well known result of Fisher and Grinstein [31 ■ This 
difference is not yet clear, but the accuracy of the mea- 
surements , as well as, the interactions from the systems 
can generate this deviation. The present results confirms 
the exciting idea of the existence of a metallic state in 
the insulator-superconductor transition which is driven 
by disorder, which in fact is not favorable for the su- 
perconducting state. We expect that this intermediate 
metallic state to appear in the the d-wave superconduc- 
tors containing nonmagnetic impurities. The model from 
Ref. [l^ contains also the influence of the temperature on 
the metallic conductivity near the critical concentration 
which destroys d-waves supercoducting state. 
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II. MODEL. 

Before the presentation of the model of an interacting 
Bose system the in a random field we start with a short 
presentation of the model proposed by Dalidovich and 
Phillips which contains the physics of the problem. 
This discussion is also relevant to make clear the validity 
and advantage of the model proposed in this section. 

A. Dalidovich-Phillips model 

The Landau action for this problem can be obtained 
using the replicas to perform the average on the disor- 
dered state. The quadratic and quartic terms describe, 
in the spin glass theory the interaction between spins 
,that appears by the average procedure on the disor- 
dered state, and can be decoupled by introducing the 
auxiliary fields Qll(k,k' ,t,t') =< S;^Sj^(k',r) > and 
^J^(k, t) =< >, where the superscriptions a and b 
represent the replica indices. A finite value of ^E*^ is 
equivalent with the phase ordering in the charge of 2e 
condensate, and < 4'° >= describe the disordered 
phase. As we mentioned in introduction for the quan- 
tum spin glass the diagonal elements of (3(k, k', r, r'), in 
the limit |t — r'| oo, are the effective order parameter. 

I 



In this Hamiltonian {Ti^ — 2m = 1), <i>(x)is the bosonic 
field, ro the control parameter (for the standard case it is 
the chemical potential ), and uq > is the bare coupling 
constant. The random field /i(x) is a Gaussian random 
variable with a Fourier transform /i(k) which satisfies: 

< h{k) >=< h*(k) >= 0, < h*(k)h(k') >= (5k,k'q (2) 

where <> indicates an average over the possible configu- 
rations of the random field and will serve to introduce a 
an effective Edwards- Anderson spin-glass parameter 
denoted by g.This Hamiltonian has been used to develop 
a functional theory [ll|, described by the action : 

= + (3) 

where: 

m 

^ J2ii^o+k^-iuJn)6ab-{q/T)dLJno]Mk)Mk) 

a,b—l k 

(4) 

In this equation k = (k, Wn) the indices a, b and m are 
the replicas indices from the standard spin-glass theory 
[l3| , method used in [ll|, [l3| to calculate the free energy 



The behavior of this parameter leads to the change of the 
critical dinamic exponent z = 1 to z = 2. The free en- 
ergy given in Q contain a Gausssian and a quartic term 
as in the standard <I>*-theory, and the contribution which 
couples the charge and the glassy degrees of freedom, and 
is considered by the authors as dominant in the bicritical 
point. However, using such a form for action the appli- 
cation of the RG method seems to be very difficult, and 
even the authors considered only the Gaussian approxi- 
mation. 

We will show that this state, studied in 0] can be mod- 
elled by an action which describes the effect of a ran- 
dom field on the interacting Bose system, using the RG 
in the quantum limit. In the following we will present 
the model which gives a similar result with the Gaussian 
model containing the corrections given by the interaction 
between bosons. The effect of disorder is contained only 
in the gaussian contribution , and we will show that in 
the lowest order the gaussian term is identical with the 
expression from [3| . 



B. Random field model 

We consider a d-dimensional Bose model in a random 
field /i(x) described by the Hamiltonian: 



(1) 

I 

of this disordered system. In the replica-trick theory the 
calculations are performed taking m — > 0. The interac- 
tion contribution has the form: 

m 

a— 1 k\ k4 

(5) 

The Gaussian propagator in the limit m — > has the 
form gP: 

GaAk) - Go{k)Sa,b + f3Gl{k)q6^^,o (6) 

where l/Tand Go(k,cj„) = {ro+k'^ + 'n\\ujn\y^ .We 
mention that the damping term linear in energy has been 
introduced in the model following the physical consider- 
ations presented in and it has a importance, be- 
cause it keeps the universality class z = 2 for the model. 

In the next section we will use the results from [HI to 
write the flow equations for the model in the limit m — > 
0. We also considered the replica symmetric case, as in 
Ref.[l3|, but we do not have a term which coupes the 
bosonic field $ and the order parameter q, which in this 
model has been considered only spatial dependent. This 



H= J (i''x$^(x)[-V^+ro]|*(x)|+^ J d'ix|*(x)|* + J d'ix[h(x)*t(x) + h*(x)*(x)]. 
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may be considered as a very poor approximation, be- 
cause we used the an free propagator with the Ohmic 
dissipative term, which impUes a time dependence of the 
effective parameter q. Anyway, even this simple model, 
which is tractable by RG- theory, contains more than the 
Gaussian model, and can give us an idea of the simple 
approximation performed in [l3l] . 



III. RENORMALIZATION GROUP EQUATIONS 

The flow equations for the RG. can be obtained using 
the recursion relations [ll| (where we take b = e') as: 



dT 



dq 
'dl 



dr 

— = 2r + uFi+ uqF2 
dl 



[d- 



z)]u--[m 



2Fa 



5u^q 



(7) 

(8) 
(9) 
(10) 



In these equations we used for the number of components 
of the bosonic field n = 2, and we will take d = 2 and 
z = 2. The functions Fi,{i = 2... 5) can be calculated pl| . 
(in this case there is a difference because of the damping 
T]) and only the function Fi{T{l) = K2[exp{l/T{l))-1]-^ 
{K2 = l/(27r)) is important because contains the relevant 
parameter T{1). The Eqs.(9-10) will be solved in the low 
temperature limit and when we can take: 8-^3-1-2^4 ~ K2 
, andFs ~ K2. 

The solutions of the Eqs.(7-8) are T{1) = Te^' and 
q{l) = qe^K Using these results we solved exactly the 



Eq.(9) and the solution is: 

u(0- 



1 



K2I + I0 + {e^ 



(11) 



where/o = 4:/K2Uo. 

The general solution of the Eq.(9) has the form 



r{l) = e^'[ro + h{T) + 4q / dl' - 



o2;' 



I' + lQ + 5q{e^' -1)' 



(12) 



where Ii (T) is given by the expression: 

We will consider the case of low temperatures when Ii 
can be neglected, and q is small, so we can approximate 
Eq.(12) by : 



r{l) 



4q / dl' 

l' + lo 



(14) 



This integral can be given by the function Ei{x) defined 
by 



Ei{xy) 



1 fV-l 

-.^y I dt 



X + lnt 



(15) 



For y — 1 and big argument Ei(x) ~ /x and we get for 
Eq.(14) the expression: 



r(0 



ro + 2g 



(16) 



Next important step is to calculate the stop scaling 
parameter from the condition r{l*) = 1. This proce- 
dure applied first in [? ] gives the possibility to study 
the influence of the temperature on the critical behav- 
ior of the transport properties in the critical region of 
the z = 2 insulator-superconductor transition. In our 
case the Gaussian fixed point T = is stable, and we 
will expect a correction given by the interaction between 
fluctuations on this behavior. Indeed from Eq.(16) we 
obtain the important result: 



{l + uoqf 



(17) 



a result which will be used to calculate the conductiv- 
ity of the system considered, at T = 0. The conductivity 
will be calculated using the Kubo formula generalized 
for the replicated action (1) which contains the influ- 
ence of disorder. In the next section we will calculate 
a{iujri) following the approximations from 13] but using 
the renormalized value for q{l) which is a relevant pa- 
rameter. This parameter is strongly dependent of the 
variable I and in this way we can fix it near the critical 
point studding the competition between disorder and the 
interaction between fiuctuations in the occurrence of the 
metallic phase. Before starting this calculation, which 
will follow the approximations from Q , we mention that 
even if our model is slightly different from that devel- 
oped by Dalidovich and Phillips [7, 's'] their approxima- 
tions remain valid for our model because are related to 
the replica-trick method, and of the propagators of the 
replicated systems, which are the same for both mod- 
els. On the other hand, the existence of the dissipation 
is essential for the occurrence of the metallic state as 
was first pointed out in Q. However, if our model ap- 
pear to be more tractable analytically by RG method, the 
main approximation of our model remain the dissipative 
form of the free propagator Go(k, CL'n).Such a form can 
be easy obtained if we adopt a model similar with this 
from 15] where the free propagator contains the scatter- 
ing effects on the non-magnetic impurities, or the cou- 
pling to an Ohmic heat bath [l^. This approximation 
is well justified by the results, and a RG-treatmcnt con- 
taining more complicated form of the action is difficult to 
be controlled. The higher order corrections mentioned in 
[T] leads the terms proportionally to u^q^ which appear 
also in our calculations, but we cannot consider these 
corrections as having the same origin. In fact the ver- 
tex corrections presented later by the authors in Ref. Q 
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leads to the same qualitative behavior for the T = 
case, the difference being the power of rg from the de- 
nominator, a result similar with the calculations from 
, performed for the case of layered quasi- two dimen- 
sional superconductor. In fact the generalization of our 
model for this case is very simple problem and we do 
not expect qualitatively new results. The origin of the 
driving parameter rg is another problem which can be 
discussed because this transition is a T = transition, 
and in fact is a Quantum Phase Transition (QPT). In the 
next section we used the results from 0, [13] for the cal- 
culation of conductivity, without giving many technical 
aspects, and using the same approximations. However, 



Following the way from [7| we transform this equation in 
the simple form: 

^(^„) = M!l£! f ^klGl{k,0)[Go{k,oyGo{k,ujn)] 

„(19) 

In this point we mention that, differently from [2] we 
considered the Z-dependence of the effective order param- 
eter q for the disordered state in order to take into consid- 
eration the influence of the interaction. This is not neces- 
sary in the case of the Gaussian approximation mention 
above. Using now the simple expression for Gq and the 
Eq.(17) for I* from Eq.(19) and q{l*) = ge"*'* we get: 

a{u = 0, T ^ 0) = const^-n^±^ (20) 
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IV. CONDUCTIVITY 

The Kubo formula for a disordered system has the gen- 
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(18) 



I 

This result shows that for uq = 0, which is Gaussian 
approximation we get the same result as in Ref [3], and 
the metallic state is robust to the interaction between 
bosons. 
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